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Updating schemes in zero-temperature single-spin flip dynamics 
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In this paper we examine the role of the so called c-parallel updating schemes in relaxation from 
disordered states to the final ferromagnetic steady state. We investigate two zero-temperature single- 
spin flip dynamics on a one dimensional lattice of length L: inflow (i.e. generalized zero-temperature 
Glauber dynamics) and outflow opinion dynamics. The varying c allows us to change the updating 
scheme from random sequential updating (c = 1/L) to deterministic synchronous updating (for 
c = 1). We show how the mean relaxation times depend on c and scale with the system size L. 
Moreover, we empirically find an analytical formula for the ratio between mean relaxation times for 
inflow and outflow dynamics. Results obtained in this paper suggest that in some sense the original 
zero-temperature Glauber dynamics is a critical one among a broader class of inflow dynamics. 

PACS numbers: 05.50.+q 
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I. INTRODUCTION 

Recently, there has been renewed interest in the pos- 
sibly different physics arising from the sequential or syn- 
chronous execution of the microscopic update rule of the 
spins in disordered systems [l[ (and references therein). 
Two most used updating schemes are - random sequen- 
tial and synchronous (parallel) updating. Within syn- 
chronous updating all units of the system are updated at 
the same time. Random sequential updating means that 
only one, randomly chosen, unit is updated at each time 
step. Probably neither a completely synchronous nor a 
random sequential update is realistic for natural systems 
0. Therefore, in the previous paper Q we introduced 
c-parallel updating, in which a randomly chosen fraction 
c of spins is updated synchronously. Varying c allows 
us to move from random sequential updating (c = 1/L) 
to deterministic synchronous updating (for c — 1). An 
analogous idea was proposed very recently by Radicchi 
et al. @ to study quench from infinite temperature to 
zero temperature of Ising spin systems evolving accord- 
ing to the Metropolis algorithm, in which the flipping 
probability associated with the iih spin is defined as Q : 



p M = mm{l,exp(-AEi/k B T)} . 



(1) 



where AEi is the energy change due to the ith spin flip. 
For AEi = we obtain pm = 1 , hence zero-temperature 
limit of this algorithm can be rewritten as a special case 
of the so called generalized zero-temperature Glauber dy- 
namics in which flips for AEi > are forbidden and 
the probability of flipping ith spin in case of AEi — 
is Wo = 10. Another example of such a generalized 
dynamics (which we call in this paper inflow dynam- 
ics) is the zero-temperature limit of the original Glauber 
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dynamics Q - the corresponding flipping probability is 
Wo = 1/2. 

It occurred that for Wo = 1 the critical value of c ex- 
ists above which the system never reaches the final ferro- 
magnetic state [2] . This critical value was determined by 
measuring the evolution of active bonds (bond is active if 
it connects two sites with opposite spins). We expect that 
in the case of generalized Glauber dynamics (Wo £ [0, 1]) 
the critical value of c depends on Wq. However, deter- 
mining precise dependence between the critical updating 
scheme c and flipping probability Wo is not the subject 
of this paper and will be considered in a future study. 
In the present paper we focus on measuring the mean 
relaxation times to the ferromagnetic steady state as a 
function of c and Wo. Moreover, apart from the inflow 
dynamics, we investigate the so called outflow dynamics, 
which was originally proposed to study opinion dynamics 
0. 

The paper is organized as follows. In the next section 
we present the motivation for this study. In the follow- 
ing sections we recall definitions of studied dynamics, as 
well as, the idea of c-parallel updating and present sim- 
ulation results. We conclude the paper by summarizing 
the results. 



II. MOTIVATION 

As Kenrick et al. Q report, about 20 years ago, for no 
good reason the customers of a local bank in Singapore 
began drawing out their money in a frenzy. The run on 
this respected bank remained a mystery until much later, 
when researchers interviewing participants discovered its 
peculiar cause: an unexpected bus strike had created an 
abnormally large crowd waiting at the bus stop in front 
of the bank that day. Mistaking the gathering for a crush 
of customers poised to withdraw their funds from a fail- 
ing bank, passersby panicked and got in line to withdraw 
their deposits, which led more passersby to do the same. 
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Soon illusion had become reality and, shortly after open- 
ing its doors, the bank was forced to close to avoid ruin. 
This is only one of many examples of social validation 
described in social psychology textbooks [1, [9|, [l(| • 

Most people feel that behaviors become more valid 
when many others are performing them. Although the 
tendency to follow the lead of our peers can drive to mis- 
guided behavior, most of the time it sends us in right 
directions, toward correct choices. On the other hand, 
it happens that people often follow others even if they 
believe that the group may be in error. 

The research of Solomon Asch in 1955 demonstrated 
that, when faced with a strong group consensus, people 
often conform even if they believe that the group may 
be in error [Tlj |. However, even a single visible dissenter 
from the group's position emboldens others to resist con- 
formity. A number of later experiments showed that an 
individual who breaks the unanimity principle reduces 
social pressure of the group dramatically [10j . This ob- 
servation was recently expressed in a simple one dimen- 
sional USDF ('United we Stand, Divided we Fall') model 
of opinion formation Q . The model was later renamed 
'the Sznajd model' by Stauffer and generalized on a 
two dimensional square lattice. 

The crucial difference between the Sznajd model and 
other Ising-type models of opinion dynamics [ID, [l4 J151 
[iH UtI Il8j ] or zero-temperature Glauber dynamics gjis 
that information flows outward from the center nodes to 
the surrounding neighborhood (so called outflow dynam- 
ics 0) and not the other way around. In spite of this 
difference, it has been suggested [H, [2(| , that outflow dy- 
namics is equivalent with zero-temperature Glauber dy- 
namics (inflow dynamics) for the Ising ferromagnet with 
the nearest neighbors interactions. 

The updating of an extended group of spins in the 
Ising model with zero-temperature Glauber kinetics m 
was considered by Newman and Stei n |2I| , Lipowski |22j 
and later on by [3, i, i, [H, H, HI, [Hll (and refer- 
enced therein). Outflow d yna mics in two dimensions was 
studied by Stauffer et alJllj2§j2iJMLand later on by 
others (for reviews see [3l[ l32i l33l. 134 l35j|. Qualitative 
differences between outflow and inflow dynamics in two 
dimensions have been shown. Under outflow dynamics 
the system always relaxes to the ferromagnetic steady 
state. On the contrary, an Ising ferromagnet has a large 
number of metastable states with respect to the Glauber 
spin-flip dynamics. Therefore, at zero temperature the 
system could get stuck forever in one of these states [24| . 
On the other hand, since there are no metastable states in 
one dimension, the only possible final states are all spins 
up or all spins down, analogously to outflow dynamics. 

Let us recall two differences between settings of inflow 
and outflow dynamics: 

f. Under zero-temperature Glauber (inflow) dynam- 
ics conformity is caused by majority, in contrast to 
outflow dynamics where unanimity is required. Of 
course, in one dimension, when a pair of spins is 
influencing others, majority is equivalent to una- 



nimity - this difference is important only in higher 
dimensions. 

2. Under inflow dynamics the center spin is influenced 
by its nearest neighbors, while under outflow dy- 
namics information flows from the center spins to 
the neighborhood. 

Nevertheless, under random sequential updating, both 
dynamics seem to be equivalent in one dimension, except 
for one value of parameter Wq = 0, which describes the 
probability of flipping the spin in lack of majority. This 
particular value Wq = corresponds to the constrained 
zero-temperature Glauber dynamics where the only pos- 
sible moves are flips of isolated spins. The system, there- 
fore, eventually reaches a blocked configuration, where 
there is no isolated spin 

In the previous paper we have decided to examine dif- 
ferences between both dynamics in one dimension under 
several types of updating [3J] . This issue seems to be quite 
important in the field of sociophysics, since several social 
experiments showed that updating plays an important 
role in opinion dynamics. For example, two groups of 
two people influence an individual stronger to conform 
than one group of four people [l^, On the other 

hand, updating schemes may play a prominent role in 
one-dimensional Ising model at zero temperature [l|, 0] • 

Probably neither a completely synchronous nor a ran- 
dom sequential update is realistic for natural systems 0] • 
Therefore, in Q we have introduced c-parallel updating 
(a randomly chosen fraction c of spins is updated syn- 
chronously) and shown that outflow dynamics is much 
more influenced by the type of updating than inflow dy- 
namics. The results in our previous paper were intrigu- 
ing, but preliminary. We looked at the mean relaxation 
time as a function of the initial fraction of randomly dis- 
tributed up-spins only for several values of c and Wq. 
Based on our previous simulations we noticed that: 

1. Within -^-parallel updating (i.e. random sequential 
updating) the relaxation is much slower for inflow 
than for outflow dynamics. 

2. Within 1-parallel updating (i.e. deterministic syn- 
chronous updating) the relaxation under outflow is 
slower than under inflow dynamics. 

3. In general, the relaxation times decay with Wq. 

4. The dependence on c is much stronger for outflow 
dynamics. For inflow dynamics the mean relax- 
ation time is almost the same for all values of c. 

As we show in this paper, the above results are valid 
only for some values of Wq and c. Generally, the pa- 
rameter space of our model is 3-dimensional and con- 
sists of: c (fraction of spins updated synchronously), Wo 
(probability of flipping a spin in lack of majority) and 
p (initial fraction of up-spins). However, in this paper 
we focus on random initial conditions with equal num- 
ber of up and down spins, i.e. p = 0.5. This value of 
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p is the most frequently chosen in the computer simula- 
tions since it corresponds to high temperature situation 
in which magnetization is equal to 0. On the other hand, 
condition p — 0.5 corresponds to the social situation in 
which public opinion is equal to zero, i.e. a democratic 
stalemate in which the society is not able to make any 
common decision. From such an initial condition, un- 
der one-dimensional zero-temperature outflow or inflow 
dynamics with random sequential updating, the system 
will eventually reach with equal probability one of two 
possible consensus states - all spins up or all spins down. 
Obviously, the mean relaxation time for p — 0.5 is larger 
than for any other value of p. For all these reasons the 
initial state with p = 0.5 is the most interesting case. On 
the other hand, focusing on this particular value of p al- 
lows us to limit the parameter space (p, c, Wo) to (c, Wo). 
In particular, we show in this paper: 

• how the mean relaxation time r scales with the 
system size L, 

• what is the exact dependence between the mean 
relaxation time and parameters c and Wo, 

• what is the ratio between the mean relaxation time 
for inflow and for outflow dynamics as the function 
of parameters c and Wq. 



III. DYNAMICS AND UPDATING 

Let us begin with recalling the definitions of inflow and 
outflow dynamics, as well as, the idea of the so called c- 
parallel updating. The system consists of L Ising spins 
Si = ±1 (i = 1,...,L) placed on the one-dimensional 
lattice with the periodic boundary conditions. We con- 
sider a quench from the infinite temperature to the zero 
temperature, and let the system then evolve using one of 
two dynamics - the inflow and the outflow dynamics. 

Generalized zero-temperature Glauber (inflow) dy- 
namics can be defined without using the idea of the en- 
ergy 3: 

Stft+l) = 

/ Si+i(t) if Si_i(t)Si+i(f) = 1, (2) 

\ -Si(t) withprob W if Si-i(t)S i+1 (t) = -1 

Similarly we can define the outflow dynamics. Recently 
[37l | , we have introduced slight modifications with respect 
to the original outflow rule: choose pair of neighbors and 
if they both are in the same state, then adjust one (in- 
stead of two) of its neighbors (chosen randomly on left or 
right with equal probability 1/2) to the common state. 
Because this way at most one spin is flipped in one step, 
while in original formulation two can be flipped simul- 
taneously, the time must be rescaled by factor i. We 
measure the time so that the speed of all processes re- 
mains constant when L — > 00, so normally one update 
takes time Here, instead, we consider also the factor 



|, so single update takes time At = tjj. Our modifica- 
tion eliminates some correlations due to simultaneous flip 
of spins at distance 3. However, if we look at later stages 
of the evolution , where typically the domains are larger 
than 2, simultaneous flips occurs very rarely. Therefore, 
we do not expect any substantial difference. Indeed, com- 
puter simulations confirmed our expectations - only time 
has to be rescaled. On the other hand, the modification 
simplifies the analytical treatment and allowed us to find 
exact formula for the exit probability [37|. Here we go 
even further and adjust left (instead of random) neighbor 
of chosen pair: 

Si(t + 1) = 

/ $+i(t) ifS i+ i(t)S i+2 (t) = l, (3) 

\-Si{t) withprob W if S i+1 (t)S i+ 2(t) =-1 

We have checked by the computer Monte Carlo simula- 
tions that such a modification does not change results in 
the case of random sequential updating. Analogously to 
the case with one random neighbor [37j , only the time has 
to be rescaled by factor 2 in respect to the original rule 
in which two spins are changed in each elementary time 
step (see Fig. [T]). We introduce such a modification to 
simplify the case with synchronous updating. Moreover, 
it should be noticed that the case of synchronous updat- 
ing for the zero-temperature Glauber dynamics is fully 
deterministic. By introducing modification described in 
Eq. ([3]) we avoid randomness in the case of synchronous 
updating for the outflow dynamics. 

Seemingly there is only one small difference between 
both dynamics - under inflow (see Eq. ©) spin Si is 
flipped according to its two nearest neighbors and 
Si+i, and under outflow (see Eq. Q) according to its 
neighboring pair S%-\-i,Si+2. This small difference has 
consequences even for random sequential updating - for 
Wo = we have constrained inflow dynamics where the 
only possible moves are flips of isolated spins and the 
system, therefore, eventually reaches a blocked configu- 
ration, where there is no isolated spin Q. On the other 
hand, under outflow dynamics with Wo = the system 
will always eventually reach consensus or, in other words, 
ferromagnetic state. However, apart from the limit value 
of the flipping probability Wo = both dynamics seem 
to be qualitatively equivalent under random sequential 
updating. To see more differences between dynamics we 
have introduced c-parallel updating [3]. Within this up- 
dating a randomly chosen fraction c of spins is updated 
synchronously. Of course, c = 1 corresponds to deter- 
ministic synchronous updating and c = 1/L to random 
sequential updating. 

IV. SIMULATION RESULTS 

In the previous paper Q we did not check how results 
scales with the system size L. It is known that the mean 
relaxation time in one-dimensional system scales with the 
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FIG. 1: The mean relaxation times from random initial state 
consisting of p up-spins for the modified (Inn) and original 
(2nn) outflow dynamics in one dimension. In the modified 
version we adjust only left neighbor of chosen pair, so at most 
one spin is flipped in one elementary step while in original 
formulation two can be flipped simultaneously. Therefore in 
case of modified version the time was rescaled by factor |. 
It should be noticed that in computer simulations time is 
measured in Monte Carlo Steps (MCS). As usual, one MCS 
consist of L elementary updating. Since we investigate the 
relaxation process we simulate the system as long as it reach 
the final state with all spins up or down. The average number 
of MCS needed to reach the final step depends on the initial 
concentration of up-spins and scales with the system size as 
L 2 analogously to the voter model [3Sl . |39| . l4p| ]. The results 
presented on the plot are averaged over 10 4 samples. 



lattice size as < r >^ L 2 for the voter model [^um|40| 
as well as for the inflow and the outflow dynamics with 
random sequential updating [3|. We decided to examine 
the scaling laws for c-synchronous updating with c > 
1/L. It occurs that in the case of inflow dynamics the 
mean relaxation time scales with the system size as < 
t >= aL 2 , where a = a(Wo) and for a given value of 
Wo is c- independent (see Fig. [5]). This results is valid 
for Wo < 0.5. However, it has been shown recently that 
for Wo = 1 there is a critical value of c = c, above which 
the relaxation time is infinite - the system never reaches 
the ferromagnetic steady state Q. Probably in general 
c = c(Wn), but this is not the subject of this paper. 
For Wo < 0.5 the system always eventually reaches the 
ferromagnetic steady state (i.e. c = for Wo < 0.5) and 
the scaling law < r >= a(Wo)L 2 is valid. 

In the case of the outflow dynamics the mean relax- 
ation time scales with the lattice size as < t >~ aL 2 
and a = a(c, Wo) (see Fig. [3]). Results presented in Figs. 
[2] and [3] confirm observation from our previous paper [3[ 
- for Wo < 0.5 outflow dynamics is more influenced by 
the type of updating than inflow dynamics. It should 
be noticed that Wo, with respect to social applications, 
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FIG. 2: The mean relaxation time from disordered to the 
final ferromagnetic steady state for one-dimensional inflow 
dynamics with c-paralel updating for c = !(*), c = 0.5(x) 
and c = I/L(o) as a function of lattice size L. The mean 
relaxation time scales with the lattice size as < r >~ L 2 . The 
results presented on the plot are averaged over I0 3 samples. 



can be interpreted as the probability of opinion change 
in case of uncertainty (lack of unanimity) and therefore 
it should be relatively small due to the observation that 
an individual who breaks the unanimity principle reduces 
social pressure of the group dramatically [l(| ■ 

Let us now examine the dependence between the mean 
relaxation time < r > and the type of updating (charac- 
terized by c) for different values of flipping probabilities 
Wo- Several interesting phenomena are seen in Figs. [4] 
and [5] In the case of the inflow dynamics (see Fig. Q|: 

1. The mean relaxation time < r > is c-independent 
for flipping probability Wo < 0.5. 

2. For Wo > 0.5 the mean relaxation time increases 
rapidly with c above certain value c — c* (e.g. for 
W = 0.6, c* w 0.7, for W = 0.8, c* » 0.45). It has 
been found recently 0] that for Wo = 1 the sys- 
tem never reaches its ferromagnetic ground state 
for c > c = 0.41, (i.e. for W = 1 and c > 0.41 
< t >= oo). We expect that c* is correlated or 
even equivalent with critical c. However, determin- 
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FIG. 3: The mean relaxation time from disordered to the 
final ferromagnetic steady state for one-dimensional outflow 
dynamics with c-parallel updating for c = 1(*), c = 0.5(x) 
and c = l/L(o) as a function of lattice size L. The mean 
relaxation time scales with the lattice size as < r >~ L 2 . The 
results presented on the plot are averaged over 10 3 samples. 



ing precise dependence between the critical updat- 
ing scheme c and flipping probability Wo is not the 
subject of this paper and will be considered in a 
future study. 

3. Generally, the mean relaxation time decreases with 
increasing Wq- However, for Wq > 0.5 this state- 
ment is valid only for c < c*. 

In the case of the outflow dynamics (see Fig. [5]): 

1. The mean relaxation time < t > depends on c for 
all values of Wq. More precisely, the mean relax- 
ation time increases with c, which is desired result 
in modeling social systems - it was shown in series 
of social experiments that two groups of two people 
influence an individual stronger than one group of 
four people [l(| . 

2. For Wq < 0.5 the mean relaxation time decreases 
with increasing Wq for whole range of c, analo- 
gously to inflow dynamics. 
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FIG. 4: The mean relaxation time from disordered to the 
final ferromagnetic steady state for one-dimensional inflow 
dynamics with c-parallel updating as a function of c. For 
Wo < 0.5 the mean relaxation time is c-independent. For 
Wo > 0.5 the mean relaxation time increases rapidly with c 
above certain value c = c* . The results presented on the plot 
are averaged over 10 3 samples for L = 100. 
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FIG. 5: The mean relaxation time from disordered to the 
final ferromagnetic steady state for one-dimensional outflow 
dynamics with c-parallel updating as a function of c. In 
case this case the dependence on c is visible for all values of 
Wq. The results presented on the plot are averaged over 10 3 
samples for L = 100. 



Let us now present dependence between the mean re- 
laxation time and flipping probability Wo (see Figs. [6] 
and . In the previous paper it was argued that in gen- 
eral, the relaxation times decay with Wo 3- However, 
this statement is valid only for c — ► For larger val- 
ues of c dependence < t(Wq) > is non monotonic for 
both (inflow and outflow) dynamics. The mean relax- 
ation time decays with Wq for Wq < Wq (c) and grows 
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above this value. This is very intriguing result but, due to 
partially synchronous updating scheme, we were not able 
to provide any analytical treatment that helps to under- 
stand this non-monotonic behavior. As we have already 
mentioned, we expect the critical value of c*(Wq), and 
analogously Wq(c), above which the relaxation time is 
infinite. However, due to infinite relaxation (simulation) 
time, the critical values should be determined by mea- 
suring the evolution of active bonds instead of relaxation 
time. This will be studied in a future. 
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FIG. 6: The mean relaxation time from disordered to the 
final ferromagnetic steady state for one-dimensional inflow 
dynamics with c-parallel updating as a function of the flip- 
ping probability Wo- The results presented on the plot are 
averaged over 10 3 samples for L — 100. 
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The main aim of this paper is to compare the mean 
relaxation times for the inflow and the outflow dynamics. 
Therefore we have decided to determine the ratio between 
mean relaxation times for inflow and outflow dynamics. 
For Wo < 0.5 this ratio be can be approximated by linear 
function of updating scheme c (see Fig. [SJ; 



< T > in 

< T > out 



(4) 



where dependence between a and Wq can be approxi- 
mated from simulation data (see Fig. [TU)) as: 



aW + 1 



0.43. 



Finally: 



< T >out 

thus < t >i n =< t > out for: 
_ aWo + 0.5 
aW + 1 



(aW + 1) c + (aW + 1.5) . 



{a w 0.43) 



(5) 
(6) 

(7) 



which in case of Glauber dynamics, i.e. Wq =0.5 gives 
c w 0.48. 
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FIG. 8: The ratio between mean relaxation times for inflow 
and outflow dynamics as a function of c for Wo < 0.5. In this 
case the ratio Ti„/T ut decays linearly with c. 



For Wo > 0.5 ratio < r >i n / < r > ou t is non- 
monotonic (see Fig. [5]). It decays with growing c for 
c < c*(Wo) and above grows rapidly with c. Interest- 
ingly this special value of flipping probability Wo = 0.5 
corresponds to original Glauber dynamics 0]. 



FIG. 7: The mean relaxation time from disordered to the fi- 
nal ferromagnetic steady state for one-dimensional outflow 
dynamics with c-parallel updating as a function of the flip- 
ping probability Wo. The results presented on the plot are 
averaged over 10 3 samples for L = 100. 



V. SUMMARY 

In this paper we examine the role of the so called o 
parallel updating schemes in relaxation times from disor- 
dered states to the final ferromagnetic steady state. We in- 
vestigate two zero-temperature single-spin flip dynamics 
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FIG. 9: The ratio between mean relaxation times for inflow 
and outflow dynamics as a function of c for Wo > 0.5. No 
linear dependence is seen in this case. 




FIG. 10: For Wo < 0.5 the ratio Tin/Tout decays linearly with 
c and it can be approximated as Ti n /r ut = —o,c + (a + |), 
where a ~ Wq 1 . Solid line represents power-law fit a = 
0.43WJ" 1 + 1. The same plot in log-log scale (inset) shows 
clearly the power-law dependence. 



on a one dimensional lattice of length L: inflow (i.e. gen- 
eralized zero-temperature Glauber dynamics) and out- 
flow opinion dynamics. Under random sequential updat- 
ing, both dynamics seem to be equivalent in one dimen- 
sion, except for one value of parameter Wq = 0, which 
describes the probability of flipping the spin in lack of 
majority. Therefore, we have decided to examine dif- 
ferences between both dynamics in one dimension under 
generalized c-parallel updating. This issue seems to be 
particularly important in the field of sociophysics. In Q 
we have introduced c-parallel updating (a randomly cho- 



sen fraction c of spins is updated synchronously). An 
analogous idea was proposed very recently by Radicchi 
et al. [2| to study quench from infinite temperature to 
zero temperature of Ising spin systems evolving accord- 
ing to the Metropolis algorithm, in which the flipping 
probability associated with the ith spin is defined as [H : 



p M = mm{l,exp(-AEi/k B T)} , 



(8) 



where AEi is the energy change due to the ith spin flip. 
For AEi = we obtain pu = 1, hence zero-temperature 
limit of this algorithm can be rewritten as a special case 
of the so called generalized zero-temperature Glauber dy- 
namics defined by Eq.((2]), with Wq — 1. It occurred that 
for Wo = 1 the critical value of c exists above which 
the system never reaches the final ferromagnetic state 
Q ■ This critical value was determined by measuring the 
evolution of active bonds. We expect that in the case of 
generalized Glauber dynamics (Wo S [0,1]) the critical 
value of c depends on Wq. However, determining pre- 
cise dependence between the critical updating scheme c 
and flipping probability Wq was not the subject of this 
paper and will be considered in a future study. In the 
present paper we focused on measuring the mean relax- 
ation times < t > to the ferromagnetic steady state as a 
function of c and Wq for inflow and outflow dynamics. 

We have found the critical value of Wo = 0.5 (which 
corresponds to the original Glauber dynamics 0), below 
which the system eventually reaches the final ferromag- 
netic steady state with probability 1 for any value of c. 
Moreover, for Wq < 0.5: 

1. the mean relaxation time scales with the system 
size as power law < t >~ L 2 for both dynamics, 

2. the mean relaxation time decreases with increasing 
Wq for both dynamics, 

3. the ratio between the mean relaxation times for 
inflow and outflow dynamics can be approximated 
by the linear function: 



< T >out 



= - (aWo + 1) c + (aW a + 1.5) . 



(9) 



4. the mean relaxation time < r > is c-independent 
for the inflow dynamics in contrast to the outflow 
dynamics. 

For Wq > 0.5, in the case of the inflow dynamics, the 
critical value of c = c* exists. Above this value the prob- 
ability of reaching the final ferromagnetic state decreases 
to zero. Simultaneously, the mean relaxation time in- 
creases rapidly with c above the critical value c — c* . 

Let us stress here, that for social applications results 
for Wq < 0.5 are particularly important, since Wq can be 
interpreted as the probability of opinion change in lack of 
majority. A number of social experiments showed that an 
individual who breaks the unanimity principle reduces so- 
cial pressure of the group dramatically [icij and therefore 
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we assume that Wq is relatively small for social appli- 
cations. Results obtained in this paper suggest that the 
outflow dynamics is more adequate for modeling social 
phenomena than the inflow dynamics, since for Wq < 0.5 
the outflow dynamics is c-dependent in contrast to the 
inflow dynamics. 

On the other hand, inflow dynamics is much more ad- 
equate for physical applications. Moreover, it should 
be noticed that value Wq = 0.5 corresponds to origi- 
nal Glauber dynamics @. Therefore, results obtained in 
this paper suggest that in some sense the original zero- 
temperature Glauber dynamics is a critical one among 
a broader class of inflow dynamics. The existence of 



the critical values of c and Wq in the c-parallel updat- 
ing scheme for the generalized zero-temperature Glauber 
dynamics are very intriguing results and deserve deeper 
analysis. 
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